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Abstract 

We present results for light hadrons composed of both degenerate and non-degenerate 
quarks in quenched lattice QCD. We calculate masses and decay constants using 60 
gauge configurations with an 0(a)-improved fermion action at (3 = 6.2. Using the 
p mass to set the scale, we find hadron masses within two to three standard de- 
viations of the experimental values (given in parentheses): mx* = 868 + g MeV 

(892 MeV), = 970 MeV (1020 MeV), m N = 820 ^ MeV (938 MeV), 

m A = 1300 MeV (1232 MeV) and m n = 1650 +™ MeV (1672 MeV). Direct 
comparison with experiment for decay constants is obscured by uncertainty in current 
renormalisations. However, for ratios of decay constants we obtain fx/ fx = 1-20 + ^ 

(1.22) and f^/f p = 1.13 (1.22). 
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1 Introduction 



Within the quenched approximation, it is currently possible to study lattice QCD numerically 
in a box of linear size around 2 fm, with a lattice spacing of less than 0.1 fm, corresponding to 
a cutoff above 2 GeV. Although extrapolation to the chiral limit is still a necessary ingredient 
when u and d quarks are involved, such lattices allow the direct simulation of hadrons 
containing s quarks, at the cost of fixing one additional mass parameter. A significantly 
wider range of physical quantities thereby becomes calculable, with which to probe the 
reliability of lattice QCD. 

In this paper, we extend an earlier study [2j of light hadron masses and decay constants 
to include the effects of S"?7(3)-flavour-symmetry breaking. The earlier study was based 
on 18 configurations, and the correlation functions were evaluated for hadrons composed 
of degenerate quarks, for five different values of the quark mass, using both the Wilson 
and the clover fermion actions. The results presented in this paper were obtained from 
our complete data set of 60 configurations, using the clover fermion action, for three of the 
five previously-used quark masses. In order to study flavour-symmetry-breaking effects, we 
construct mesons using all possible quark-mass combinations. 

In the Wilson formulation, the bare quark mass, m, is given in terms of the hopping param- 
eter k by 



K cr it is the value of the hopping parameter at zero quark mass, which is taken to be the point 
at which the mass of the pseudoscalar meson, mp, vanishes and the quark and antiquark are 
degenerate. With each quark flavour, we need to associate a value of k corresponding to its 
experimentally observed mass. It is a good approximation to take the physical light quark 
(u and d) masses to be zero, i.e. to set n u = Kd = « C rit- Here we are interested in computing 
the effect of S , f/(3)-flavour-symmetry breaking on the spectrum and decay constants, and 
so we need to associate a non-zero mass with the strange quark. The corresponding k s 
can be determined, for example, by first extrapolating the vector meson mass, my(/«i, K2), 
to K\ = K2 = Kcrit and then fitting the data for the ratio mp(K\, K2)/"T.y( /? crit ) ^crit) to 
some function of the two quark masses, extrapolating in K\ to n CTit and using « 2 to fix the 
ratio to the experimental value of m 2 K /m? p . Alternatively, k s may be determined from the 
degenerate-quark data alone, avoiding the second chiral extrapolation, by using the ratio 
ma/m p . 

This procedure requires some assumption about how the hadron masses depend on the quark 
masses. In most previous calculations (see, for example, [@, 0, f|, ||, |6|), it has been assumed 




(1) 



that the pseudoscalar meson mass obeys the PCAC relation 

/J_ J 1 

\2ki 2k 2 K crit 



mp(«i, k 2 ) = 6p ( h — ), (2) 



and that the vector meson mass obeys 

mv{Ki, K2) = a v + by ( — + — —) . (3) 

The original work of Martinelli et al. [0] supported this assumption, albeit on the basis of 
rather limited statistics, by combining quark propagators computed for different values of k. 
More recent studies of strange hadrons [§, [| have similarly utilised non-degenerate quarks, 
and this is the procedure that we adopt here. 

We fit our data for mp, my, fp, etc., to the following function of the two quark masses, mi 
and 777-2, 

^2 / \ ^3 11 / \ 

ai + — (777 2 + mi) + — \m 2 -mi\ (4) 

to test the assumption that a 3 = in Equations (fj) and (Q). The expression (|]), or its 
equivalent with = 0, defines a plane through the data, as shown in Figure IB. By means 
of a combination of extrapolation and interpolation using such fits, we are able to calculate 
the masses and decay constants of the K, K* and <fi mesons and the mass of the Q baryon, 
in addition to those of the usual light hadrons. 




Figure 1: Schematic of the fit planes defined by Equation (f|). 

In Section 2, we summarise our statistics and present details of our fitting criteria and proce- 
dures. Our main results are given in Section 3, where we describe the chiral extrapolations, 



the interpolation to the strange quark mass, and the masses and decay constants in physical 
units. Section 4 contains our conclusions. 



2 Computational Details 

2.1 Statistics 

We have analysed 60 gauge field configurations at an inverse coupling (3 = 6.2 on a lattice 
of size 24 3 x 48. The gauge configurations and quark propagators were produced on the 
64-node i860 Meiko Computing Surface at the University of Edinburgh. The SU (3) gauge 
fields were generated using the Hybrid Over- Relaxed algorithm, defined in reference . The 
gauge configurations are separated by 2400 sweeps, beginning at configuration 16800. The 



quark propagators were calculated using an 0(a)-improved clover action [10. I I 



Sp = Sp 7 - E q{x)F lliy (x)cr llu q(x). (5) 
Sp 7 is the standard Wilson lattice action, 
Sf = E{?(z)^) - K E [?(*)(! - l^(x)q(x + q(x + + ^)Ul(x)q(x)} \ (6) 

and F^v is a lattice definition of the field strength tensor. We have computed propagators 
at three values of k, 0.14144, 0.14226 and 0.14262, using an over-relaxed minimal residual 
algorithm with red-black preconditioning and point sources and sinks. Although there are 
advantages to using smeared sources and/or sinks to extract ground-state properties from 
2-point functions constraints imposed by other parts of the UKQCD programme did 



not allow us this option. 

We construct correlators for mesons composed of quarks of flavours 1 and 2 using the fol- 
lowing local interpolating fields: 

P = <?l75<?2 (7) 

A A = gi7475<?2 (8) 

Vi = qi%q 2 , (9) 

and correlators for baryons composed of degenerate quarks using 

N = e abc (u a C l5 d b )u c (10) 

A M = e afec ( M a C 7 yK. (11) 



For the vector meson, we average our correlators over the three polarisation states, for 
the nucleon we average the 11 and 22 spinor indices of the correlator, and for the A we 



project out the spin-| component and average over the four spin projections. Our quark 
propagators incorporate the rotations required to ensure that using these interpolating fields 
yields 0(a)-improved correlators. This and our computational procedure are described in 
detail in reference 0. 

Except where explicitly stated otherwise, the errors quoted in this paper are purely statistical 
and are calculated according to the prescription: 

• create 1000 bootstrap samples from the original dataset of 60 configurations by ran- 
domly choosing, with replacement, 60 configurations per sample; 

• for each bootstrap sample, perform all the mass fits and extrapolations as for the 
original data; 

• obtain the errors on a given quantity from the 68% confidence limits of the correspond- 
ing bootstrap distribution. 



2.2 Fitting Procedure 

We construct 2-point meson correlation functions from quark propagators with all com- 
binations of the three k values. We perform least-% 2 fits to the zero-momentum, time- 
symmetrised time-slice correlators to single cosh functions. For the pseudoscalar channel, 
we fit over the time range t — 14 to 22 for all k combinations. For the vector channel, we 
use the fitting range t = 13 to 23 for all except the heaviest degenerate-K case, where we use 
t = 15 to 23. For the baryons, we construct 2-point correlation functions only for the degen- 
erate cases. We fit the appropriate average of the forwards and backwards, zero-momentum 
time-slice correlators to single exponential functions, choosing the time ranges t = 16 to 22 
for the nucleon and t — 16 to 21 for the A. We take account of time correlations in the 
least-x 2 fits. 

We carried out an extensive investigation of the most appropriate fitting ranges, before 
arriving at the above choices. With reference to the effective mass plots in Figures |2] and |^, 
we fixed t max to be as large as possible and reduced t m i n until the x 2 /dof showed a significant 
increase. In this way, we attempted to fit as many time slices as possible. We followed the 
reasoning of reference [jnj in fitting our data as far out as possible to avoid contamination 



from excited states at earlier times. 

In Figure ^ we show the variation of the degenerate-quark vector meson mass estimates 
and of the corresponding x 2 /dof with the position of the first time slice, t m i n , in a variable 
window, fixing t max = 23. We observe a slow monotonic decrease from t min = 14 onwards, 
at all k values, although there is no significant variation in x 2 /dof from t min = 13 and the 
mass estimates agree within statistical errors for the two lightest k values. The case of 
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Figure 2: Effective mass plots for the pseudoscalar, vector, nucleoli and A at k — 0.14144. 
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Figure 3: Effective mass plots for the pseudoscalar, vector, nucleoli and A at k — 0.14262. 
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Figure 4: Fit stability plots for the vector meson. The top graph shows the variation of 
fit mass with t m i n , fixing t max = 23. The bottom graph shows the corresponding x 2 /dof. 
The different symbols refer to different k values: squares = 0.14144, diamonds = 0.14226, 
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m v 


statistical 


systematic 


0.14144 


0.395 


0.376 


0.389 


+0.007- 0.006 


+0.006- 0.013 


0.14226 


0.351 


0.312 


0.343 


+0.009 - 0.007 


+0.008- 0.031 


0.14262 


0.335 


0.284 


0.319 


+0.014-0.013 


+0.016-0.035 



Table 1: Highest and lowest fit masses for the vector meson with degenerate quarks, based 
on an analysis of the fit regions 13-23 to 18-23. We take the difference between the highest 
(lowest) mass and the best-fit mass as a measure of the systematic error. 











error 
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high 

m N 


m N 


m N 


statistical 


systematic 


0.14144 


0.573 


0.568 


0.573 


+0.015-0.007 


+0.000- 0.005 


0.14226 


0.475 


0.452 


0.462 


+0.020 - 0.014 


+0.013-0.010 


0.14262 


0.392 


0.372 


0.372 


+0.024- 0.016 


+0.020- 0.000 



Table 2: Highest and lowest fit masses for the nucleon, based on an analysis of the fit regions 
16-22 to 19-22 for the heaviest and 14-22 to 19-22 for the others. We take the difference 
between the highest (lowest) mass and the best-fit mass as a measure of the systematic error. 

the heaviest k value is more problematic given the small statistical errors. We adopt the 
compromise choice of t min = 15 for this, which is consistent with a wider range of other 
choices of t m ; n . We attempt to quantify in Table [l] the systematic error arising from the 
choice of fit range, by taking the lowest and highest mass estimates from fits with values 
of £ min acceptable according to the above criteria. We note that the monotonic decrease 
in our mass estimates with increasing t m i n is reflected in the asymmetry of the systematic 
error estimates towards lower mass values. However, beyond this, we feel unable to quote 
quantitative estimates of this error in our results. 

In Figure |5| and Table ^ we show the corresponding data for the nucleon. These provide 
convincing evidence of plateaux at all three k values, beginning at t m i n = 16 for the heaviest 
and at t m i n = 14 for the two lightest masses. We note that the large upper error bars in 
the latter data at t m i n = 14 and 15 may indicate the sensitivity of the bootstrap sampling 
to the tail of the excited-state contributions. For this reason, we take t min = 16 for all three 
k values. The systematic error estimates in Table are well within the statistical errors, 
giving us confidence in our choice of fit range. 
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Figure 5: Fit stability plots for the nucleoli. The top graph shows the variation of fit mass 
with t m in, fixing t max = 22. The bottom graph shows the corresponding x 2 /dof. The different 
symbols refer to different k values: squares = 0.14144, diamonds = 0.14226, circles = 0.14262. 



We calculate the pseudoscalar decay constant from the ratio 
£aj<Ai(aj,*)pt(o)) f p7Tlp 



tanhmp(L 4 /2 — t), 



(12) 



Zx(P(x,t)Pt(0)) Z A (0\P\P) 
fitting to the same time slice range as in the fit to the pseudoscalar correlator. We use the 
parameters from this latter fit to fix mp and the matrix element (0\P\P). We find that this 
particular ratio of correlators gives the cleanest signal from which to extract fp fl2|]. We 
determine /y by fitting to 



3m 



v 



££<^0M)^(o)> 

j=l X ZZj VJV 



e -m v L t /2 coshmy ( Lt /2-t). 



(13) 



Here, Za and Zy are the factors required to ensure that the lattice currents obey the correct 



current algebra in the continuum limit [14, 15 



3 Results 

3.1 Masses and Decay Constants in Lattice Units 

In Table |3| we present the masses and decay constants of the vector and pseudoscalar mesons 
calculated from both degenerate-quark and non-degenerate-quark correlators. The x 2 /dof 
are all satisfactory, being generally between 0.5 and 2. The degenerate-quark data may be 
compared with our results from the first 18 configurations, presented in 0, where our fit 
range was 12 - 16 for all hadrons. Our new estimates are within la of the 18 configuration 
estimates. The errors are reduced by a factor of approximately 2. 

As noted previously 0, @, experimental data suggests that the hyperfine splitting, m v — mp, 
should be only weakly dependent on the quark masses for light hadrons, and this we observed 
within large errors. Our higher-statistics results quoted in Table [J[ are shown in Figure [f| 
The new data for both degenerate and non-degenerate quarks is entirely consistent with our 
earlier results, but with significantly smaller errors. 

In Table |] we present the masses of the nucleon and A computed using degenerate quarks 
only. Again, the x 2 / dof of the fits is acceptable. Our estimates for the nucleon mass are about 
2a below our estimates based on the first 18 configurations [0, with only slightly smaller 
errors. It is clear from Figure [| that t m i n = 12 is not in the plateau region when using 60 
configurations. We have used the increase in sample size primarily to reduce the systematic 
error in the mass due to contamination by excited states. However, this may not be the only 
effect. We observe increasing fluctuations in the nucleon mass with decreasing quark mass, 
and fluctuations which decrease the nucleon mass tend to dominate the statistical average. 
The influence of such fluctuations in the present analysis may be more pronounced because 
we fit the nucleon correlator further from the source. 
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Table 3: Masses and decay constants in lattice units of mesons composed of degenerate and 
non-degenerate quarks. 
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Figure 6: Vector-pseudoscalar meson mass splitting for both degenerate- and non-degenerate- 
quark data. 
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Table 4: Masses in lattice units of baryons composed of degenerate quarks. 



K 


mp/mv 


m N /mv 


0.14144 


0.77 +J 


1.47 +\ 


0.14226 


0.62 + _\ 


1.35 +; 


0.14262 


0.52 + _\ 


1.17 + " 

— b 



Table 5: Mass ratios from degenerate-quark data, used in the Edinburgh plot. 



Apart from at the highest quark mass, where our new estimate is 3a lower, our estimates 
for mA agree with our previous results. The errors have not decreased, probably because of 
the extended fitting range. 

We have looked for evidence of correlations between successive configurations, by varying 
the bin size in a jackknife error analysis. This did not reveal any significant effects in the 
hadron time slice correlators. 

In Figure ^ we show the Edinburgh plot for our degenerate-quark data. The corresponding 
mass ratios are given in Table ^. The apparently alarming fall of the data points may not 
indicate any discrepancy with experiment, as we shall see that the chirally-extrapolated 
value of m^jmy is 1.07 * , within 2a of the experimental value, as may be deduced from 



Table [10]. Indeed, our higher-statistics results are only 1 — 2a below our earlier result using 
local sources and sinks Jl]. We attribute this difference to our new, lower estimates for the 
nucleon mass. We remark that were the systematic errors in the vector meson mass, discussed 
above, to be included somehow, the tendency would be for the upper and right-hand error 
bars to increase. 
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Figure 7: Edinburgh plot for degenerate-quark data. 
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Table 6: Fit parameters for masses, decay constants and ratios, using the fit form described 
in equation (f|), (A) with a 3 unconstrained, and (B) with a 3 = 0. 



3.2 Quark-Mass Dependences 

We fit our lattice estimates for the meson masses, decay constants and their ratios, for all 
six quark-mass combinations, to the expression (H) with a 3 , the coefficient of \rri2 — 

(A) unconstrained and (B) constrained to be zero. In Table ^| we present values for the fit 
coefficients. For all the quantities, the values obtained for a 3 in fit (A) are consistent with 
zero to within one or two standard deviations. The x 2 /dof for all the fits are satisfactory, 
although mostly a little larger for fit (B), and in each case the two fits give completely 
consistent values for a% and This is the numerical evidence that these physical quantities 
depend only on the sum of the quark masses. As there is no theoretical justification for 
a 3 7^ 0, nor any support for this from our data, from hereon we use only the results of fit 

(B) . 



3.3 Chiral Extrapolations 



Firstly, we fit our estimates of the pseudoscalar meson mass for all six quark-mass combina- 
tions to the form in Equation (|2|), and obtain K cr \ t from extrapolating the fit in both k's to 



m 2 p(K crit , Kent) = 0. This gives 

n ciit = 0.14315 +1 (14) 

in good agreement with, although significantly more accurate than, our estimate from 18 
configurations |§. In Figure || we present the plot of m 2 P versus l/2n e s from which we derive 
K cr it, defining an effective n as 

1 - ( 1 + M (15) 
K e ff \2ki 2k 2 / 

The fact that both the degenerate- and non-degenerate-quark data agree well with the fit 
is graphical evidence of our claim in the previous section that the pseudoscalar meson mass 
depends only on the sum of the quark masses, and confirms the observation of reference @ . 
The chiral extrapolation of my to K CTit using Equation ([$]) is also shown in Figure [5]; again, 
we conclude from the good agreement between the data and the fit that m v depends only 
on the sum of the quark masses. 

We present, in Figure |9], the chiral extrapolations of the nucleon and A. It is evident that 
for the nucleon the quality of the linear fit is rather poor. This is supported by the fact that 
the x 2 /dof for this fit is 5.1, compared to the value 2.4 obtained for the pseudoscalar meson 
fit and 0.63 obtained for the vector meson fit. Negative curvature of the nucleon has been 
observed before (see for example || |16[). A linear fit to our data for m 2 N gives a smaller x 2 /dof 
of 2.7, as well as a substantially smaller estimate for the nucleon mass in the chiral limit. 
However, because there is no theoretical justification for this choice of extrapolation, and 
because we cannot reliably compare different choices with only three data points, we quote 
results only from the linear extrapolation of tun, noting that there is significant uncertainty 
in this procedure. We see no such problem with the chiral extrapolation of m^, although 
the errors are larger. 

The inverse lattice spacing, a -1 , in physical units, obtained from each of the p, nucleon and 
A masses, is given in Table |7|. Throughout the remainder of this paper, we will use m p 
to set the scale. In so far as there is good agreement between the scales from the string 
tension and m p , our results in physical units are not especially dependent on the chiral 
extrapolation. The corresponding values for the nucleon and A masses in physical units are 
given in Table [It]. 



3.4 Determination of the Strange Quark Mass 

We calculate the value of k s from the fit to the ratio m 2 P (Ki, K2)/m 2 p for all six quark-mass 
combinations, extrapolating K\ to K cr ; t , and using k<i to match the ratio to its experimental 
value, m 2 K /m 2 = 0.413, giving 

K2 = Ks = 0.1419 + _\. (16) 
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;ure 8: Chiral extrapolations of the pseudoscalar and vector meson masses. 
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;ure 9: Chiral extrapolations of m N and m& for degenerate-quark data only. 
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Table 7: Chir ally-extrapolated lattice masses and corresponding scales compared with the 
scale from the string tension 0. 

This value agrees well with that obtained from degenerate-quark data with smeared sinks in 
reference [0. We note that m s lies between two of our quark masses, so that our strange- 
quark results are obtained by interpolation, an intrinsically more robust procedure than 
extrapolation. 

Our value of k s corresponds, in the lattice regularisation, to a strange quark mass in physical 
units of 

m s (a) = Ya(^-^ L )= 82 - 8 s MeV < ^ 
where we have used m p to determine the lattice spacing in physical units. From this we can 
determine the renormalised strange quark mass in the MS scheme at a reference renormali- 
sation scale of \i = 2 GeV: 

m^(2GeV) = 109 +JJ, (18) 

where we have used the perturbative value of the renormalisation constant relating m s (a) 
and (//) Ul7| , together with the effective coupling defined in reference ||18|| . This result 
agrees with the conclusion from previous lattice evaluations of the strange quark mass (see 
reference [O for a simulation with Wilson fermions and reference [EDI for one with the 



clover action at (3 = 6.0), that m^ s (2 GeV) is at the lower end of expectations [21], 22|. The 
value given in Equation ( |T8| ) can be compared to 100 ± 6 MeV quoted in reference [I9| and 
89 ± 9 MeV quoted in reference pUJ . 



3.5 Spectrum and Decay Constants in Physical Units 

We present the results for meson masses and decay constants, extrapolated/interpolated to 
the physical k values, in Table |8|, using the parameters of fit (B). As indicated in Table ||, 
the ratios of the decay constants need to be multiplied by renormalisation constants before 
they can be compared with experimental numbers. The perturbative estimates for these 
renormalisation constants |15|, using the effective coupling, are 

Z A ~ 0.97, Z v ~ 0.83. (19) 



rriK* 



ni 



m 



K 



Itt/Za 

Ik/Za 

fx/ fn 



y(f P zv) 

Wk*Z v ) 

U/ fp 

fx* I fp 



Ul (m p Z A ) 
f K / (m p Z A ) 

fx/ ( m K* Z a) 



0.251 
0.326 
0.364 
0.075 



+ 13 
-11 
+ 13 
-12 
+ 12 
-11 
+ 6 
-5 



0.040 
0.048 
1.20 



f i 

-2 
f 1 
-2 
-3 
-2 



0.380 
0.359 
0.337 
1.13 
1.06 



+ 9 
-16 
+ 6 
-11 
+ 4 
-7 
+ 2 
-3 
+ 1 
-2 



0.142 
0.165 
0.148 



+ 6 
-9 
+ 7 
-8 
+ 4 
-6 



Table 8: Extrapolated/interpolated values of meson masses and decay constants in lattice 
units, calculated using the fit (B) parameters in Table ||. 



Incorporating these values for the renormalisation constants and using the lattice scale from 
m p in Table [7], we obtain the meson masses and decay constants in physical units presented 
in Table |[ 

The results for m K *, m K « — m K and m^, although 2 — 3a below experiment, provide sup- 
port for our determination of k s . However, they are open to the interpretation that the 
vector meson masses may be slightly underestimated relative to the pseudoscalar meson 
masses, as a result of the suppression of spin-splittings in the quenched approximation. Our 
result for txik* is significantly more precise than that obtained with the standard Wilson 
action at (3 = 6.0 by Lipps et al. § who quote 930(40) MeV, and by Loft & DeGrand § 
who obtain 761(122) MeV. It is comparable with the value 896(17) MeV obtained using a 
renormalisation-group-improved action by Iwasaki p. Lipps et al. and Iwasaki both use m$ 
to determine k s , based on the assumption that the is pure ss, and so we are only able to 
compare our result for with that of Loft & DeGrand, who give 868(114) MeV. Loft & 
DeGrand determine n s from the T] s , which is the hypothetical pure ss pseudoscalar meson, 





lattice estimates 


experiment 




670 + *°MeV 

— 1U 


"686 MeV" 


rriK* 


868 + ;?MeV 

— o 


892 MeV 




970 ^°MeV 


1020 MeV 


m 2 K * — m 2 K 


0.53 + J (GeV) 2 


0.55 (GeV) 2 


u 


102 + ^MeV 


132 MeV 


f K 


123 + ^MeV 

— 6 


160 MeV 




0.316 + J 3 


0.28 


1/fK* 


0.298 + l 






0.280 + l 

— 


0.23 


U/m p 


0.138 +g 


0.172 


fK/m p 


0.160 + l 

— o 


0.208 


f K /m K * 


0.144 + i 

— 


0.179 



Table 9: Values of meson masses and decay constants in physical units, using the scale from 



whose theoretically-expected mass of 686 MeV [ffl agrees well with our calculation. 

Martinelli and Maiani |5j noted that the ratio (m^* — m p ) / {m 2 K — m 2 ) typically gives a 
value for the inverse lattice spacing which is lower than that obtained using other physical 
quantities. The ratio is estimated by assuming that my and mp are linear in the sum of the 
quark masses, Equations (0) and (|3]), so that: 

my(«i, k 2 ) = a v + ^-m P (Ki, k 2 ) (20) 
bp 



and 

m K , - m p _ m v (K s , K crit ) - my(/t crit , Kent) _ by_ 
m\-m\ mp(K s , K cri t ) - m 2 P ( K CTit , K cri t ) b P ' 



(21) 



We obtain 2.3 * ^ GeV, in good agreement with our low-statistics results 0. This is con- 
sistent with the scales we obtain from other physical quantities, although it remains on the 
low side. 

Our results for the pseudoscalar decay constants are 5 — 7cr below their experimental values. 
This has been noted recently in simulations using the standard Wilson action ||23|| . These 



authors argue that a smaller value of the decay constant is to be expected in the quenched 
approximation than in the full theory, as a consequence of the smaller wavefunction at the 
origin. The discrepancy may also be partly due to our use of the perturbative value for Z^. 
Recent non-perturbative calculations at (3 = 6.0, and at a single value of the quark mass, put 
the value of Z& ~ 1.09(3) ||24|| , somewhat higher than the perturbative value. Thus, it may 
be hoped that a full non-perturbative evaluation of Za at (3 = 6.2 will raise our estimates of 
the decay constants, bringing them closer to the experimental values. The uncertainty in the 
renormalisation constant is removed in the ratio fx/ fit = 1-20 ^> which agrees well with 
the experimental value of 1.22. This suggests that, although we are working in the quenched 
approximation, we obtain correctly the dependence on the strange quark mass. We cannot 
attribute the problem with the overall normalisation to quenching until we have a precise 
non-perturbative determination of the axial current renormalisation. 

Experimentally, the ratio of the pseudoscalar decay constant to the vector meson mass is 
fairly insensitive to the S[/(3)-navour-symmetry breaking. In Figure |TD| we plot this ratio 
versus mp for both the degenerate- and non-degenerate-quark data. The slope agrees well 
with that of the experimental data and we note that an increase in Za of order 25% would 
give excellent agreement between the two. 

Our estimates for the vector meson decay constants, 1/ f p and 1//,/,, given in Table |9], lie above 
the experimental values. This is shown for all our data in Figure [D]. The discrepancy is small 
for the p, but it becomes significant for the 0, although the slope of the data is consistent with 
experiment. The sign of the discrepancy is opposite to that expected from the suppression 
of the wavefunction at the origin by the quenched approximation, suggesting that other 
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Figure 10: fp/my against mp in physical units using the scale from m p . Experimental 
points for the pion and kaon are also shown. 



effects are important. The difference is less likely than in the case of the pseudoscalar decay 
constant to be due to our use of the perturbative value of the renormalisation constant, 
because at least at (3 = 6.0 Zy is known to be close to its non-perturbative value [23J. Thus, 
we suspect that significant discretisation and/or finite- volume errors are present. 
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Figure 11: 1/fv against my in physical units using the scale from m p . Experimental points 
for the p and mesons are also shown. 



For the baryon data in Table |], we perform extrapolations and interpolations in the single 
k value to obtain our predictions for the physical masses. The A and Q baryons both have 
J p 



| , so by interpolating our A fits to k s we can obtain ijiq, as shown in Table [L0[ We 



are encouraged by the good correspondence between our computed values and the experi- 
mental masses presented in the table. However, the low confidence we have in the linear 
chiral extrapolation of the nucleon means that we cannot attach much significance to the 
comparison between the computed and experimental values for tjin- The extrapolation of 
ira does not have this problem, and the interpolation to k s gives a value of very close to 
the experimental number. This indicates that our calculation of k s can be applied sensibly 
to the baryon sector. Our value for itlq is to be compared with that of Lipps et al. [|| who 
quote 1650(150) MeV, Loft & DeGrand § who quote 1512(144) MeV, and Iwasaki [§ who 
quotes 1737(77) MeV. 





extrapolated values 


experiment 




lattice units 


MeV 




m N 


0.31 +;| 


820 «° 


938 MeV 


m A 


0.50 ti 


1300 ™ 


1232 MeV 


m n 


0.62 + _\ 


1650 ™ 


1672 MeV 


m n /m N 


2.0 


f 2 
-2 


1.78 



Table 10: Physical values of the baryon masses using m p to set the scale, including the ratio 
m n /m N . 



Our determination of the spin-| baryon mass provides us with an alternative means of 
determining k s , by interpolating our lattice data to the physical value of m^. Using m p 
to set the scale, we obtain in this way k s = 0.1417 * , ra r)s = 740 ^ MeV and = 

I 2Q 

1010 _Z MeV. Although there remains an implicit dependence of these estimates on the 
chiral extrapolation, through our use of m p , this can be avoided by, for example, taking the 
scale from the string tension. It is evident from Table [7| that this would give similar values. 
Thus, this method of determining k s depends only on the mild assumption that the baryon 
mass varies smoothly with quark mass close to the strange quark mass. The agreement, 
within the somewhat larger statistical errors, with the value of k s obtained from rriK, in 
Equation (|i~6"D , and with the corresponding mass estimates in Table [|, provides a further 
check on our determination of n s . Finally, we note that assuming that the <ft is pure ss does 
not in practice permit the determination of k s from hadrons composed solely of s quarks, 
because the numerical data for the ratio m^jmy^K, k) is only weakly dependent on k. 



4 Conclusions 

Despite working in the quenched approximation, our simulations using the clover action 
at (3 — 6.2 demonstrate good agreement for meson masses and decay constants with the 
dependence on light quark masses, including S'[/(3)-flavour-symmetry breaking, expected 
from chiral perturbation theory. In particular, we find evidence from simulations with non- 
degenerate quarks that the dependences on the bare quark masses, 

mp(«i,«2) = h p{^- + lT — ) ( 22 ) 

m v (K 1 ,K 2 ) = av + bv(-^— + 7T ~~) > ( 23 ) 



hold at least for quark masses up to that of the strange quark. Our results for the strange- 
particle spectrum are encouragingly close to the experimental values, and the dependence 
of both the pseudoscalar and vector decay constants on m s is consistent with experiment. 
The only serious problem in the meson sector shows up in the actual values of the decay 
constants. For the pseudoscalar case, the sign of the discrepancy is consistent with the 
expectation that the wavefunction at the origin is reduced in the quenched approximation, 
but whether the magnitude is entirely attributable to this effect or is due in part to our use 
of the perturbative value for the axial current renormalisation is not known. Our results 
for the vector decay constants are more difficult to interpret, as the discrepancy is in the 
opposite direction, which may signal the presence of other effects. It is clearly important to 
calculate the current renormalisations non-perturbatively at (5 — 6.2 for at least two values 
of the quark mass. 

We have not explored S'f/(3)-flavour-symmetry breaking in the baryon sector, but only 
present results for degenerate quarks. Compared with the chiral extrapolations for the 
other hadrons, a linear chiral extrapolation for the nucleon is not well supported by our 
data, and leads to a nucleon mass which is lower than is typical of quenched simulations to 
date, although in better agreement with experiment! On the other hand, our results for the 
A and Q are less problematic. They are also in good agreement with experiment, and the 
result for the Q, along with our results for mesons, encourages us to believe that we have 
good control of the strange-quark physics that can be obtained from 2-point functions. 
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